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A REVIEW OF THE UNKNOWN INPUT OBSERVER WITH EXAMPLES 


SAM NAZARI 


Abstract. We consider the Unknown Input Observer and focus on three examples found in the literature. 


1. INTRODUCTION 

This memo considers the design of observers for a class of linear dynamic systems in which system uncertainty can 
be modeled as an additive unknown disturbance term in the dynamic equation: 

(l.I) x{t) = Ax{t) + Bu{t) + Ed{t) 

y{t) = Cx{t) 

where x{t) G M" is the state vector, y(t) G is the output vector, u(t) G K’' is the known input vector and d(t) G 
is the unknown input or disturbance vector. If the unknown input distribution matrix, E, is not full column rank, 
then the rank decomposition: 


Ed{t) = EiE2d(t) 

can be applied where, Ei, is a full column rank matrix and E 2 d{t) can be considered as a new unknown input. 
The term, Ed(t), may be used to model additive disturbances or model uncertainties that are not known a priori. 
Typical uncertainties are noise, plant nonlinearities, model reduction errors, parametric uncertainties and difficult 
to model interconnection terms in large scale systems. While not considered in this memo, a disturbance term may 
also appear in the output equation: 

y{t) = Cx{t) + Eyd{t) 

Additionally, this memo does not consider the Du{t) term in the system output equations. For the sake of brevity 
and without loss of generality, we will not treat systems with this term. Note that the development of the theoretical 
material in this memo closely follows Chapter 3 of [1]. 

2. Unknown Input Observer Theory 

The problem of designing observers for a linear system with both known and unknown inputs has been studied 
since the 1970’s [D, m, m- The motivation for studying this problem is that in practice, for a variety of reasons, 
many plants are modeled with disturbance terms as shown in Eq o When modeled in such a way, traditional 
observers with the Luenberger structure become difficult to implement since they use all of the input signals in order 
to compute an estimate of the state vector. This limitation renders them ineffective for many practical applications. 
Contrastingly, the observer structure considered in this memo assumes no a priori knowledge about the unknown 
inputs, rendering it more amenable for use in fault detection applications. Many investigators have addressed this 
problem with varying success and we refer the reader to Section 3.2 of [T] for a broad literature review. The class of 
observers examined in this memo is based on the Unknown Input Observer (UIO) scheme proposed by Chen, Patton 
and Zhang. 

Definition 2.1. An observer is defined as an unknown input observer nil] for the system described by Eq O if 
its state estimation error vector, e{t), defined as: 

(2.1) e{t) = x{t) — x{t) 

approaches zero asymptotically regardless of the presence of unknown inputs, d(t), in the system. Furthermore, the 
structure for a full order UIO is given by the dynamic system: 

z{t) = Fz{t) + TBu{t) + Ky{t) 
x{t) = z{t) + Hy{t) 
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unknown input 



Unknown Input Observer (UIO) 


Figure 1. UIO Structure where K = Ki + K 2 . 


where x G M" is the state estimate, z € K" is the state of the full-order dynamic observer, and F, T, K, H are matrices 
to be designed for achieving unknown input decoupling. A block diagram of the UIO is shown in Figure [l] From 
the block diagram, one can see that the UIO is essentially a dynamic system that decouples the state estimation 
dynamics from the disturbance term in the original system. By expanding e(t) one obtains: 


e{t) = {A- HCA - KiC)e{t) + [F-{A- HCA - KxC)] z{t) 
+ [K 2 -{A-HCA-KiC)]y{t) 

+ \T-{I- HC)]Bu{t) + {HC - I)Ed{t) 

and it is easy to see that in order to make the estimation error a function of Fe{t): 


(2.3) e{t) = Fe{t) 

the equations: 


(2.4) 

0= {HC -I)E 

(2.5) 

T = I- HC 

(2.6) 

F = A- HCA- KiC 

(2.7) 

K2 = FH 


must hold. If all the eigenvalues in F are stable, then e{t) will approach zero asymptotically. Note that Eq 2.3 


IS 


not a function of E or d. Therefore, the estimation error approaches zero independently of the disturbance terms, 
achieving the desired decoupling of the state estimate from the unknown disturbance inputs. 


Theorem 2.2. The necessary and sufficient conditions for the system in Eq \2.2\ to be an UIO for the system defined 
in II. II are: 

(1) rank{CE) = rank{E) 

(2) (C, Ai) is a detectable pair, where Ai = A — E[{CEffiCE]~^{CEffiCA 


Proof. Please see [I] 


□ 


According to [4], the disturbance estimate can now be obtained, if needed, by the relation: 


( 2 . 8 ) 


d = {CE)^[^ - CAx - CBu] 


Next we provide a step by step design procedure for the UIO that will be utilized in the examples that ensue. 
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2.1. Design Procedure for UIOs. We will now describe a procedure that can lead to effective UIO design. It can 
be shown that Eqn 2.4 is solvable if and only if the condition rank{CE) = rank{E) is met [T]. Therefore, the rank 
condition rank{CE) = rank{E) must first be checked. If this condition is not met, then a UIO does not exist. If 
the condition is met, then we may compute the following observer matrices: 


H = E[{CEfCE]-\CEf 
T = I-HC 


A, = TA 


Next, the observability of the pair (C,Ai) must be checked. As shown in Eq 2.3 an important step in designing a 
UIO is to stabilize F — Ai —KiC. If the pair {C, Ai) is detectable, then this can be achieved by using pole placement 
to choose the appropriate gain matrix Ki. It can be shown that the observability of the pair {C,Ai) is equivalent 
to the observability of the pair (C, A) [5]. Therefore, if the pair (C, Ai) is observable or at least detectable then the 
gain, Ki, in Eq |2.6| can be obtained for the UIO. Hence, if (C, Ai) is observable then a UIO exists and Ki should be 
computed using pole placement. 


On the other hand, if the pair (C, Ai) is not observable, then one must construct a transformation matrix, P, for 
the observable canonical decomposition of the system: 


PAiP-i 


All 0 

Ai 2 A22 


CP-^ = [C* 0] 


All e 
C* G 


where ni is the rank of the observability pair (C, Ai), and (C*, An) is observable. If any of the eigenvalues of A 22 
are unstable, then a UIO does not exist. If the eigenvalues of A 22 are stable, then (C, Ai) is detectable and we may 
select ni desirable eigenvalues and assign them to An — K^C* using pole placement to obtain K^. Next compute 


K,=p-^K^ = p-\Klf 

where K'^ can be any (n — nl) x m matrix, because it does not affect the eigenvalues of F. Finally, compute the 
observer matrices: 

P = Ai - KiC 
K = Ki+K2=Ki + EH. 

In order to clarify the design procedure we consider some numerical examples next. 


3. Numerical Examples 

3.1. Third order system with no inputs. The first example considered is taken from page 76 of [T]. There, the 
following parameter matrices are used: 


(3.1) 


A = 


-I 1 
-I 0 
0 -I 


First we check that rank{CE) = rank{E) 

rank(C*E) 

raiik(E) 


0 

0 

-I 

= I: 


■ 




, c = 

■ 1 

0 

0 ■ 

0 

0 

1 






E=[ -1 0 0 ] 


cuis = ans = 


1 


1 


The rank is equal to I as needed. The observer matrices are computed next: 
H = E*inv((C*E)’*(C*E))*(C*E)’ 

T = eye(3)-H*C 
AI = T*A 

H = T = AI = 


1 0 0 0 0 
0 0 0 1 0 
0 0 0 0 1 


0 0 0 
-10 0 
0 -1 -1 
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Next, the observability of the pair {Ai,C) must be checked. If it is full rank {rank{obsv{Ai,C)) = 3), then we may 
apply pole placement: 
rank(obsv(Al,C)) 

cins = 

3 

Now we proceed to apply pole placement in order to place the observer poles at [—2, —10, —5]. 

K1 = place(A>,C’,[-2,-10,-5]) ’ 

K1 = 

6.7231 -3.7816 

7.5581 -11.3458 

-3.0543 8.2769 

Finally, the observer F and K matrices are computed: 

F = A1-K1*C 
K = K1 + F*H 

F = K = 

-6.7231 0 3.7816 0 -3.7816 

-8.5581 0 11.3458 -1.0 -11.3458 

3.0543 -1.0000 -9.2769 0 8.2769 



Figure 2. Simulink block diagram of the dynamics subsystem. 
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Figure 3. Simulink block diagram of the UIO. 


Let us examine the Simulation results with a disturbance occurring at t = 10 seconds. This is showin in Figure 
The initial conditions of the dynamic system as: 


xl_0 = 100; 
x2_0 = -100; 
x3_0 = 1; 

In Figure]^ one is able to see the state trajectories in the solid lines. The dotted ’x’ represent the estimate from the 
UIO. For the given disturbance and the chosen observer pole locations, the observer seems to converge (remove all 
estimation error) within five seconds. Figures]^ andshow the block diagram in detail. Note that we have chosen 
not to estimate the disturbance itself. All the files for this simulation can be found in the 

C:\Users\sqn4594\Documents\FDIR 


folder. Specifically, the file 


bookEx.m 


and the Simulink model 


exOne.xsl 


contain everything. 
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Dynamic system state and UlO estimate 



Figure 4. State response and UIO estimate for example 


3.2. Car Suspension System. The next example is a quarter car model. The system matrices are: 


r 0 

1 

0 

0 1 


r 0 

0 1 


r 1 

0 

0 

0 1 


\ ^ 1 

— ks 

— bs 

ks 

bs 


0 

1000 


0 

1 

0 

0 


1 

mb 

mb 

mb 

mb 

B = 

mb 

c = 

E = 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

1 

ks 

bs 

-ks-kt 

— bs 


kt 

-1000 


0 

0 

0 

1 


1 

L mw 

mw 

mw 

mw J 


L mw 

mw J 







where uib, represents the car body and the mass, rriw, represents the wheel assembly. The spring ,ks, and damper, 
bs, represent the passive spring and shock absorber placed between the car body and the wheel assembly. The spring, 
kt, models the compressibility of the pneumatic tire. The variables Xb-, x^ and r are the car body travel, the wheel 
travel, and the road disturbance, respectively. The force, fs, which is applied between the body and the wheel 
assembly, is controlled by feedback. We assume that the disturbance force may enter at any of the state variables, 
not just the road disturbance. In MATLAB, we assign the following values to the parameters: 

7o-System Physical Constants-7 

mb = 300; 7 kg 

mw = 60; 7 kg 

bs = 1000; 7 N/m/s 

ks = 16000 ; 7 N/m 

kt = 190000; 7 N/m 

Therefore, the system matrices become: 

7 -System Matrices (It is a SIMO system)- 7 

A = [ 0 1 0 0; [-ks -bs ks bs]/mb ; ... 

0001; [ks bs -ks-kt -bs]/mw]; 

B = [00; 0 10000/mb ; 0 0; [kt -10000]/mw]; 

C = eye(4); 

E=[l;l;l;l]; 
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Figure 5. Free Body Diagram of the Quarter Car Suspension model (idealization). 


A = 





B = 


l.Oe+03 

* 




l.Oe+03 * 


0 


0.0010 

0 

0 

0 

0 

-0.0533 


-0.0033 

0.0533 

0.0033 

0 

0.0333 

0 


0 

0 

0.0010 

0 

0 

0.2667 


0.0167 

-3.4333 

-0.0167 

3.1667 

-0.1667 

C = 





E = 


1 

0 

0 

0 


1 


0 

1 

0 

0 


1 


0 

0 

1 

0 


1 


0 

0 

0 

1 


1 



Now we can obtain an optimal gain for control of the suspension system by using LQR below. We will not dwell too 
much on this part of the design work since our main focus is the UIO design. 


7,-LQR design-7, 

Q=[.25 0 0 0;0 4 0 0;0 0 1 0;0 004]; 

R=50*eye(2); 

[K1 1 s]=lqr(A,B,Q,R); 

In order to design the UIO, we first check to see if the rank{CE) = rank{E) = 1 


rank(C*E) 

raiik(E) 


ans = ans = 


1 


1 
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Figure 6. Root leve Simulink Block diagram architecture of the car suspension system in example two. 


Next we compute the observer matrices 


H = E*inv((C*E)’*(C*E))*(C*E)’ 
T = eye(4)-H*C 
A1 = T*A 

f=zeros(4,4) 
v=[-4 -4 -4 -4] 

F=diag(v) 

kl=inv(C)*(A-F-H*C*A) 

k2=F*H 

k=zeros 

k=kl+k2 



A 


Figure 7. Quarter car suspension system dynamics implementation in Simulink. 
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Figure 8. Simulink implementation of the Unknown Input Observer for the quarter car suspension model. 
H = T = 


0.2500 

0.2500 

0.2500 

0.2500 

0.7500 

-0.2500 

-0.2500 

-0.2500 

0.2500 

0.2500 

0.2500 

0.2500 

-0.2500 

0.7500 

-0.2500 

-0.2500 

0.2500 

0.2500 

0.2500 

0.2500 

-0.2500 

-0.2500 

0.7500 

-0.2500 

0.2500 

0.2500 

0.2500 

0.2500 

-0.2500 

-0.2500 

-0.2500 

0.7500 


l.Oe+03 * 


-0.0533 

-0.0026 

0.8450 

0 

-0.1067 

-0.0069 

0.8983 

0 

-0.0533 

-0.0036 

0.8450 

0 

0.2133 

0.0131 

-2.5883 

-0 


0031 0 0 0 0 

0064 0 0 0 0 

0041 0 0 0 0 

0136 0 0 0 0 



F = eye(-4) 


kl = 

l.Oe+03 * 

-0.0493 -0.0026 0.8450 

-0.1067 -0.0029 0.8983 

-0.0533 -0.0036 0.8490 

0.2133 0.0131 -2.5883 


k2 = 4X4 matrix ’- 1 ’ 


0.0031 
0.0064 
0.0041 
-0.0096 


k = 
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Estimation Error for Car Suspension System with Disturbance 



Figure 9. UIO results for the car suspension model. 


l.Oe+03 * 




-0.0503 

-0.0036 

0.8440 

0.0021 

-0.1077 

-0.0039 

0.8973 

0.0054 

-0.0543 

-0.0046 

0.8480 

0.0031 

0.2123 

0.0121 

-2.5893 

-0.0106 


Next we form A1 and check observability 
raiik(obsv(Al ,C) ) 

ans = 


4 


Since the observability matrix has full rank, we can go on and simulate the system to obtain the dynamic estimator 
performance. The initial conditions are set to be: 

X_0 = [-1;10;3;5] 

X_0 = 


-1 

10 

3 

5 


Figure shows the state estimation error for the UIO used in the car suspension model. The observer seems to 
converge within 1.5 seconds. Figures shows the implementation of the UIO in Simulink. All the files for this 
simulation can be found in the FDIR folder. Specifically, the file 

carsus.m 

and the Simulink model 
car_sus.xsl 
contain everything. 
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Figure 10. The free body diagram for the two mass and spring system. 


3.3. Coupled Mass-Spring System. Consider the coupled spring mass system shown in Figure [TOl The coupled 
mass spring system is a well studied mechanical system in control theory. The input to this system is the sinusoidal 
motion of the end of the rightmost spring, and the output is the position of each mass, <71 and ( 72 - The disturbance 
input is also shown in the figure as e. In state space form, the equations of motion are: 



9i 


0 

0 

1 

0 ■ 


0 







■ 0 ■ 


qi 

, ^ = 

0 

0 

0 

1 

, B = 

0 

, c = 

■ 1 

0 

0 

0 ■ 

, E = 

0 

X = 

9l 

-2k 

m 

k 

m 

— c 

m 

0 

0 

0 

1 

0 

0 

0 


. 92 . 


1 

3k 

-2k 

m 

0 

— c 

m 


k 

. m . 







1 


In MATLAB, we assign the following values to the parameters: 
7„-System Parameters- "L 

k =1 7. 

m =1 7. kg 

c =17. 


k = 


1 1 


1 


to obtain the following state space formulation: 


7. -State Space Formulation- 

A =[0010; 


0 0 0 1 ; 

-2*k/m k/m -c/m 0; 
k/m -2*k/m 0 -c/m] 

B = [0;0;0;k/m] 

C = [1 0 0 0;0 1 0 0;0 0 1 0;0 0 0 1] 
D =0 

E = [0;0;0;1] 


A = B = 

0 0 10 0 
0 0 0 1 0 
-2 1-10 0 

1 - 20-1 1 


■7. 


C = D = E = 

1 0 0 0 0 0 
0 10 0 0 

0 0 10 0 

0 0 0 1 1 


In order to control the placement of the masses in the system, we utilized LQR. The following MATLAB code 
accomplishes this objective: 
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controller 


Figure 11 . The root level Simulink block diagram for the coupled mass spring dynamic system. 

7.-LQR Control-7. 

Q=[100 000; 

0 100 0 0 ; 

0 0 100 0 ; 

000 100] 

R=1 

[Klqr 11 ss]=lqr(A,B,Q,R) 

Q = R = Klqr = 

100 0 0 0 1 -1.3620 10.4615 2.0165 10.0419 

0 100 0 0 

00 100 0 

0 0 0 100 

Now we are ready to begin the UIO design. The first step is to check that rank{CE) = rank{E): 

7o-Step 1: rank(CE) = rank(E) ?= 1-7o 

rank(C*E) 

raiik(E) 

ans = ans = 

1 1 

Next we are ready to compute the first set of observer matrices: 

7o-Step 2: Compute observer matrices-7o 

H = E*inv((C*E)’*(C*E))*(C*E)’ 

T = eye(4)-H*C 
A1 = T*A 


0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 1 


10 0 0 
0 10 0 
0 0 10 
0 0 0 0 


0 0 10 

0 0 0 1 

-2 1-1 0 

0 0 0 0 


H = 


T 


A 1 
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Figure 12. The mass spring system dynamics implemented in linearized form in Simulink. 


Now that Ai has been obtained, we must check the observability of the pair (C, Ai)-. 

7,-Step 3:Check (C,A1) rank- °L 

rank(obsv(Al,C)) 


ans = 


4 


Since the pair (C, Ai) is full rank, it is observable. Therefore, we can use pole placement for the observer gain Ki: 

7o-Step 4: Choose observer poles-7. 

Kl = placeCA’,C’,[-2,-10,-5,-3])’ 


2 0 10 

0 10 0 1 

-2140 
1-202 


Kl 





















































14 


SAM NAZARI 



Figure 13. The mass spring system dynamics implemented in linearized form in Simulink. 


and to conclude the observer design we compute the F and K matrices: 


7o-Step 5: Finish observer design-7, 

F = Al-Kl=t=C 
K = K1 + F*H 


-2 0 0 0 

0 -10 0 0 

0 0-5 0 

-12 0-2 


2 0 10 

0 10 0 1 

-2140 
1-200 


The simulation results are considered next. The top level Simulink block diagram is shown in Figure 11 The mass 
spring system dynamics are implemented as shown in Figure |12| The UIO implementation is shown in Figure 13 
The state estimation error is shown in Figure [Td] The following initial conditions were used for the simulation. 


xl_0 

= -1 

x2_0 

= -5 

x3_0 

= 1; 

x4_0 

= 1; 
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State Estimation Error for coupled mass spring system 



Figure 14. The state estimation error for the mass spring system using the UIO. 

4. Conclusion 

This memo begins by motivating the UIO paradigm for state estimation. A main advantage of this approach is to 
decouple the plant disturbance inputs from the state estimation process. The UIO design procedure was reviewed 
and illustrated with three common examples found in control theory. The estimation error for the coupled mass 
spring system, the car suspension system and the linear system with no inputs was considered. In each example, it 
was shown that this error asymptotically approaches zero by utilization of a UIO. 
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